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ABSTRACT 

We present the discovery of a brown dwarf or possible planet at a projected separation of 1.9" = 
29 AU around the star GJ 758, placing it between the separations at which substellar companions 
are expected to form by core accretion (~5 AU) or direct gravitational collapse (typically >100AU). 
The object was detected by direct imaging of its thermal glow with Subaru/HiCIAO. At 10-40 times 
the mass of Jupiter and a temperature of 550-640 K, GJ 758 B constitutes one of the few known 
T-type companions, and the coldest ever to be imaged in thermal light around a sun-like star. Its 
orbit is likely eccentric and of a size comparable to Pluto's orbit, possibly as a result of gravitational 
scattering or outward migration. A candidate second companion is detected at 1.2" at one epoch. 
Subject headings: planetary systems — stars: low-mass, brown dwarfs — techniques: high angular 
resolution 



1. INTRODUCTION 

While the extrasolar planets and brown dwarf com- 
panions currently known from direct imaging mark sig- 
nificant discoveries, it is important to note that none 
of these systems present scenarios analogous to a Solar- 
like system. Almost all imaged substellar companions 
either fe ature extreme orbital separations (hu ndreds of 
AU, e.g. lLuhman et al.ll2007t iKalas et al.ll2008D. star-like 
rathe r than planet-like temperatures (e.g. iMarois et all 
2008), or host stars at the extreme ends of the mass 
spectrum (A- and l ate M-type, e.g. IKalas et"aLl 120081 : 
iChauvin et al.l l2005h . The closest match is probably 
Gl 229 B, a brown dw arf orbiting an Ml star at 40 AU 
(|Nakaiima et alj I1995D . This shows how a remarkable 
gap still exists between recent direct-imaging discoveries 
and the exploration of Sun-like systems. 

In this work, we present the discovery of a substellar 
companion with a mass of 10-40 Mj up , a temperature of 
550-640K, and a projected separation of 1.9" = 29 AU. 
The estimated semi-major axis range, overlapping with 
our own Solar System's outer planet orbits, along with 
the Sun-like parent star, make it a superior laboratory for 
advancing our conventional wisdom of planet or brown 
dwarf formation among solar- like systems. 

2. OBSERVATIONS AND DATA REDUCTION 

GJ 7 58 is a G9-type star, located at a distan ce of 
15.5 pc (|Grav et all 120031 : IPerrvman fc ESAIH997 "). Its 
mass a nd radius are about 0.97 M© and 0.88 i?©, respec- 
tively (|Takeda et al.ll200"7t) . Owing to its proximity and 
solar-like characteristics, it has been surveyed for plane- 
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tary companions with the ra dial velocity technique (e.g. 
iGretherfc Lincwcavcr 2006) , but no such companions 
have yet been reported. Furthermore, it has been stud- 
ied in search of infrared excess indicating the presence 
of a circumstellar debris disk, but with null results (e.g. 
IKospal et afll2009h . 

We first detected GJ 758 B with the HiCIAO high- 
contrast imaging instrument (|Hodapp et alJl200Sl ) m an- 
gular differential imaging mode on Subaru telescope on 
May 3, 2009, with a field of view of 20" and a pixel scale 
of 9.5 mas. The images were taken in pupil-stabilized 
mode in the near infrared (H band, 1.6 /xm), where sub- 
stellar obj ects are expected to be bright with th ermal 
radiation (jBaraffe et al. 2003; B urrows et al.ll2006f). The 
star is known to have nine field objects (jCarson et all 
2009). Data reduction of the 10 exposures of 15 s re- 
vealed a hitherto unknown tenth object in close separa- 
tion with 5 a confidence. A follow-up observation of 46 
exposures of 9.7s was obtained at the next commission- 
ing run on August 6, 2009, in which the object was suc- 
cessfully rediscovered with 8.5 a confidence. The time be- 
tween epochs was long enough to confirm common proper 
motion with GJ 758 and thus to establish a gravitation- 
ally bound system. Furthermore, the higher quality of 
the second observation revealed another faint signal of 
5.6cr at an even smaller separation, whose nature re- 
mains unknown until its physical existence and proper 
motions can be proven with a follow-up detection at a 
later epoch. Both runs had excellent weather conditions 
(0.5" natural seeing in H band). 

We used a loc ally optimized combina tion of images al- 
gorithm (LOCL lLafreniere et al1l2007af ) to maximize the 
efficie ncy of the angular differential imaging technique 
(ADI, IMarois eFail 120061 ). The signal-to -noise maps de- 
rived from the final images are shown in Fig. [T] 

3. RESULTS 

3.1. Proper motion analysis 

Located only 15.5 pc away from the Sun, GJ 758 B 
exhibits strong parallactic and proper motion, providing 
a powerful tool to distinguish physical companions to the 
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Fig. 1. — Signal-to-noise maps of the discovery images for 
GJ 758 B after ADI data reduction, for the May 2009 (top) and 
August 2009 (bottom) observations. The greyscale stretch is 
[— lcr, 5cr], The S/N ratio is calculated in concentric annuli around 
the star. GJ 758 B is marked with solid arrows. A possible second 
object in the August images is pointed out with dashed arrows. 
The white plus sign marks the location of the host star GJ 758; 
the black disks designate the regions in which the field rotation arc 
is insufficient for ADI. 

target star from unrelated background stars. In a time 
series of images centered on GJ 758, a companion should 
remain bound to its host while the background stars drift 
along the predicted background trajectory in unison. 

Since the detections of GJ 758 B are only separated by 
three months, the parallactic and proper motions are of 
order 0.1" (~10 HiCIAO pixels). The 7 background stars 
visible at both epochs were used to fine-tune the distor- 
tion correction of the images. As illustrated in Figure [H 
this yields a tight cluster of proper motion vectors for 
the background stars, with a standard deviation on the 
order of the pixel scale. GJ 758 B, on the other hand, 
is found to share common proper motion with its host 
star, confirming that they form a gravitationally bound 
system. The companion's observed motion deviates from 
the proper motion distribution of the background stars 
by 10 <7, thus a chance alignment can be excluded. 
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Fig. 2. — Proper motion analysis of GJ 758's field objects. The 
circles show the change in the positions relative to GJ 758 in the 
May and August images for all 8 objects visible in both images. 
One-pixel error bars are drawn. The curved arrow designates the 
predicted parallactic and proper motion for a background star of 
GJ 758 between those epochs. GJ 758 B demonstrates common 
proper motion with its host star as well as some orbital motion 
(dotted arrow), whereas the other objects conform with the behav- 
ior expected from background stars. The size of the PSF resolution 
element and the detector pixels is provided for reference. 

The precision of the distortion correction is 0.6 pix- 
els in the central region, of the same order as the esti- 
mated precision of the centroiding of 5 a sources. We 
therefore adopt a 1 pixel = 9.5 mas error bar for our 
astrometric measurements of GJ 758 B. The tightly con- 
strained spread of the 7 background stars, which suffer 
from greater residual distortion (1.3 pixels over the whole 
detector) and their own proper motion (typical value 0.5 
pixels in 3 months), proves this a valid assumption. 

GJ 758 B is measured to move by (-7.6±9.5, 24.0±9.5) 
mas along the (RA, Dec) axes relative to GJ 758 between 
the two epochs. This represents the orbital motion of 
the companion around its host. No proper motion test 
is available for the possible second companion since it is 
detected only in one epoch. 

3.2. System age 

The age of the star is an important quantity to de- 
rive characteristics for the companion, including the mass 
and temperature. However, since GJ 758 is a field main- 
sequence solar-like star with no known connection to any 
co-moving stellar association, the age is very difficult 
to con strain. The isochronal analysis in iTakeda et al.l 
(|2007f ) places the object at an age of 0.7 Gyr, with an up- 
per limit of 3.8 Gyr at 68% confidence. The super-solar 
meta llicity of GJ 758 ([Fe/H]= 0.14 (IHolmberg et all 
[20091 ) or [Fe/H]= 0.22 (IKospal et all 12009m might sup- 
port such a relatively young age compared to the Sun, 
but may also just reflect a locally metal-rich natal 
environment . However, the activity - rotati on-age cal- 
ibration in iMamaiek fc Hillenbrand! (|2008f ) yields an 
age o f 6.2 Gyr given th e rotational period of 39.0 
days (Wright et al.l I2004T) and chrom ospheric activity 
log r£„ = -5.06 (IDuncanet al.lll991l). The higher ac- 
tivity log i?HK = -4.94 given in lWright et all (|2004fl even 
results in 8.7 Gyr. The difference is due to monitoring 
during different parts of the stellar activity cycle. The 
discrepancy with respect to the isochronal dating could 
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TABLE 1 



Observed quantities of GJ 758 B and "GJ 758 C" . 




GJ 758 B 


LiJ too O 


Photometry on August 


6, 2009: 




App. H mag. 


19.26 ± 0.16 


18.47 ± 0.24 


Abs. H mag. 


18.30 ± 0.16 


17.51 ± 0.24 


Contrast (Amag) 


14.51 ± 0.16 


13.72 ± 0.24 


Contrast (10~ 6 ) 


1.57 ± 0.18 


3.24 ± O.oo 


Astrometry on May 3, 2009: 




Proj. sep. (") 


1.879 ± 0.005 




Proj. sep. (AU) 


29.12 ± 0.08 




Position angle (°) 


197.77 ± 0.15 




Astrometry on August 6, 2009: 




Proj. sep. (") 


1.858 ± 0.005 


1 i qq _i_ n nn^ 


Proj. sep. (AU) 


28.80 ± 0.08 


18.42 ± 0.08 


Position angle (°) 


198.18 ± 0.15 


219.16 ± 0.08 


Assumed age 


mass temp. 


mass temp. 


(Gyr) 


(Mjup) (K) 


(Miup) (K) 


0.7 


10.3 549 


11.7 631 


2.0 


16.6 592 


20.4 679 


4.5 


28.6 623 


35.0 715 


6.2 


34.3 624 


41.0 717 


8.7 


39.6 637 


46.5 733 



Notes. The conversions from flux to mass and effective tem per- 
ature are based on the CDND models by Baraffc ct al. (2003). 



possibly be explain ed with long-term activity cycles (e.g. 
Uanson eT al. 2008) for what concerns chromospheric ac- 
tivity, but the rotation is harder to explain in that con- 
text. The age from rotat ion alone, according to the cali- 
bration in iBarnesI (|2007f l. is 5.5 Gyr. 

Hence, the different age indicators for GJ 758 are dis- 
crepant. Taken together, for the purpose of our analy- 
sis, we conservatively chose the rotation-activity age es- 
timate with the higher activity measurement of 6.2 Gyr 
as the baseline case for physical interpretation. For the 
error bars, we considered the full non-overlapping range 
of 0.7 to 8.7 Gyr. 

3.3. Physical properties of GJ 758 B 

The photometry of GJ 758 B and the candidate second 
companion are listed in Table [TJ These values are based 
on the August 2009 data, since GJ 758 B is blended with 
a positive radial background structure in the May 2009 
data (see Fig.[T]). The mass and effective te mperature are 
derive d on the basis of the C0ND models by Baraffc ct al. 
(2003) co mmonly used in pub l ications on sub s tellar com- 
panions (jMarois et alJ [2001 IChauvin et alJ [20051 ). A 
sample of system ages is assumed, covering the range 
of 0.7-8.7 Gyr discussed in the target properties section. 

Flux measurements were performed in flat circular 
apertures of 5 pixels in diameter. The loss of flux due 
to angular smearing in the pupil-stabilized exposures or 
partial subtraction in the ADI data reduction was care- 
fully assessed and compensated. 

For the baseline age of 6.2 Gyr, GJ 758 B is found to 
have a mass of 34.3 Jupiter masses (Mj up ) and an ef- 
fective temperature of 624 K. At the age of 0.7 Gyr, the 
mass would drop to 10.3 Mj up , placing the companion in 
the planetary regime, whereas at 8.7 Gyr it would rise 
to 39.6Afj up , making the companion a low-mass brown 
dwarf. The effective temperature, on the other hand, 
is only weakly dependent on the system age, ranging 
from 549 K to 637 K. Given this effective temperature, 



we e xpect GJ 758 B to have a T9 type spectral pro- 
file ([Leggett et alJ [20091 ). This makes GJ 758 B the 
coolest companion to a sun-like star ever imaged (cf. 
750 ± 50 K for GJ 570 D orbitin g a K4 triple system 
at 1500 AU ([Burgasser et al.ll2000h: 570 ± 25 K for Wolf 
940 B orbiting an M4 star (|Burningham et al.ll2009h : 810 
± 50 K for HP 3651 B orbiting a more sun- like K0 star 
(|Luhman et alJ [20071 )). The temperature range overlaps 
with those of the latest-type f i eld dw arfs (Dcl orme et al.l 
120081: iBurningham et al.ll2008t iWarren et al.ll2007l ). mak 
ing GJ 758 B a candidate for the coldest body outside 
the Solar System ever thermally imaged. If the second 
detected object is assumed to be a real companion, its 
mass estimates go from 11.7Mj up (also in the planetary 
regime) to 46.5 Mj up and its effective temperature from 
631Kto 733K. 

All of these temperatures are low enough to allow a 
strong presence of methane in the atmosphere, which 
can be confirmed with narrow-band differential imaging 
in the near-infrared methane absorption bands. The cur- 
rent data set, taken entirely in the broadband H filter, 
cannot by itself exclude the possibility of a white dwarf 
companion in place of a brown dwarf or planet. However, 
the measured contrast and the maximum expected age of 
8.7 Gyr render this highly unlikely. In order for the hypo- 
thetical progenitor star to live through its main-sequence 
lifetime, evolve into a white dwarf, and cool down to 
the observed H-band luminosity within 8.7 Gyr, both 
the lifetime of the star and the cooling timescale would 
have to be on the short end of the scale, implying both 
a high progenitor mass (>2M ) and a high surviving 
white dwarf mass (~1.2 M ff ). two independently unlikely 
assum ptions ()Fontaine et al J 1200 1: Holb erg fc Bergeron! 
[2006] ). 

4. DISCUSSION 

Although the relative error on the measured orbital 
motion of GJ 758 B is large, the data suffices to gauge 
the parameter space of possible orbits and derive a "best 
guess" set of orbital parameters. A custom Monte Carlo 
simulation is used for this purpose. The code generates 
a large number (10 6 ) of orbital trajectories with random 
values for eccentricity e, inclination i, argument of pe- 
riastron lu, longitude of the ascending node Q, and a 
provisional semimajor axis. The distributions are pre- 
sumed to be flat, except for the inclination, where larger 
angles are favored proportionately to sin % in order to rep- 
resent their higher geometric likelihood. The ellipse is 
then scaled such as to match the midpoint between the 
two known companion positions within the error; this 
sets the final value of the semimajor axis a. Next, the 
physical velocity vector at the midpoint is calculated for 
each orbit, projected into the image plane, and compared 
to the measured orbital motion. If it fits within the error 
bars, the orbital solution is considered valid, which holds 
for about 6% of all orbits. Finally, the ensemble of valid 
orbits is weighted according to the ratio of mean and 
current orbital velocity, (w)/w bs, to represent the statis- 
tical likelihood of finding the companion in its current 
position. 

The results are presented in Table [2] and Figure [31 
Most notably, the estimated eccentricity is very high 
(ewM = 0.691) and constrained. No orbital solutions 
below e < 0.25 match the observations within the er- 
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TABLE 2 

Estimated orbital parameters of GJ 758 B. 





Weighted median 


68% likelihood 


Semimajor axis a (AU) 


54.5 


33.9-118.0 


Eccentricity e 


0.691 


0.497-0.866 


Inclination i (°) 


46.5 


24.0-67.4 


Period P (yr) 


291 


170-658 




°- 2 Inclination i (°) 

15 30 45 60 75 90 

; ^^^^m" —. m : : i ; 

o.o I ■ ■ i i i ■ i i i i ■ i . i i i i i ■ i i i i ■ i . . ■ ■ 
50 100 150 200 250 300 

Semimajor axis a (AU) 

Fig. 3. — Statistical evaluation of the orbital solutions matching 
the observations within error bars among 10 6 orbits generated by 
a Monte Carlo simulation. The eccentricity e is plotted against 
the semimajor axis a for a selection of 1000 orbits biased according 
to statistical weight. The weight is defined as the mean orbital 
velocity divided by the velocity at the observed position, (ii)/f bs- 
The orbits are color-coded to show their inclination. The weighted 
median value for (a, e) is marked with a white plus sign. 

ror bars. The range of semimajor axes (33.9-118.0 AU) 
places the companion in the brown dwarf desert. 

The orbital separation is too lar ge to be the result of in 
situ formation by core accret ion (|Mordasini et alj|2008t 
iDodson- Robinson et alJ l2009f> . Likewise, while gravita- 
tional instabili ty (IBossI I2003D may i n principle occur at 
these ra nges ( Folt in fe Klahrl 12 009. in p rep . ) , simula- 
tions by iStamatellos fc Whitwor th (2009) suggest that 
such objects either rapidly grow into full-fledged stars 
or are scattered to larger separations. This leaves the 
possibility of formation by core accretion and subse- 
que nt outward transp ort by either gravitational scatter- 
ing (jVeras et al.ll2009l ) or outward migration as part of a 
resonant pair of co mpanions sharing a common disk gap 
(|Crida et al.H2009h . 

Scattering typically leaves companions in eccentric or- 
bits, which matches the observation. On the other hand, 
should GJ 758 C turn out to be a real companion orbiting 
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inwards from GJ 758 B, its separation would be too large 
to be the scattering partner. Conversely, it would fit the 
pair-wise migration scenario particularly well, especially 
since it would be more massive than GJ 758 B, as the 
mechanism requires. At any rate, such a configuration 
would likely be unstable on long timescales, favoring a 
young system with planet-mass companions. 

Finally, GJ 758 B could already have formed during 
the star formation phase, as a result of fra gmentation of 
the molecular cloud (e.g. iBate et al.ll2002[ ). making it a 
"failed star" . However, such brown dwarfs are generally 
ejected from the stellar system due to interactions with 
other, more massive members of the forming cluster, or 
undergo further accretion to become a full-fledged star. 
The resulting scarcity of brown dwarfs in orbit around 
stars ( conservative upper limit of 4% in lLafreniere et all 
l2007bf ) is known as the brown dw arf desert. Accordingly , 
the hydrodynamic simulations of B ate fc Bonn ell (2005) 
favor the formation of close binaries (<10AU) roughly 
equal-mass components; the few wide binaries with mass 
ratios <0.25 are formed by chance combination of two 
stars ejected from the cluster at the same time. Since 
a 0.97 Mq star like GJ 758 belongs to the most massive 
members of its birth cluster, it is unlikely to undergo 
ejection and therefore does not fit this scenario. 

5. CONCLUSIONS 

We present the discovery of a substellar companion to 
the star GJ 758 at a separation of 1.9" by angular dif- 
ferential imaging in H-band with Subaru/HiCIAO. Com- 
mon proper motion with its host is demonstrated. We 
derive ranges of 10.3-39.6 Mj up for the mass and 549- 
637 K for the temperature of the companion, classifying 
it as either a methane-bearing high-mass planet or low- 
mass brown dwarf. A candidate second companion is 
detected at a separation of 1.2", which could represent a 
migration partner of GJ 758 B. 

This unique configuration of the GJ 758 system, its 
short dynamical timescale (observable motion in only 3 
months), and its demonstrated accessibility by angular 
differential imaging, make it an excellent showcase target 
for the direct investigation of the formation and evolution 
mechanisms of planets and brown dwarfs. 
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